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Abstract. The model under consideration is a classical 2D Coulomb gas 
of pointlike positive and negative unit charges, interacting via a logarithmic 
potential. In the whole stability range of temperatures, the equilibrium statistical 
mechanics of this fluid model is exactly solvable via an equivalence with the 
integrable 2D sine-Gordon field theory. The exact solution includes the bulk 
thermodynamics, special cases of the surface thermodynamics, and the large- 
distance asymptotic behavior of the two-body correlation functions. 



PACS numbers: 52.25.Kn, 61.20.Gy, 05.70.-a 

1. Introduction 

The classical (i.e. non-quantum) equilibrium statistical mechanics deals in general 
with two basic kinds of models: discrete lattice systems and continuous fluids. In one 
spatial dimension (ID), both kinds of statistical models, considered with short-range 
as well as long-range pairwise interactions among constituents, are solvable in many 
cases In 2D, there is a large family of integrable lattice systems exactly solvable 
via the Bethe-ansatz method (see reviews ||^ and [|j). On the other hand, there was 
no exactly solved fluid in more than ID. The only partial exceptions were represented 
by 2D logarithmic models of Coulomb fluids, the one-component plasma [Q and the 
symmetric two-component plasma (or Coulomb gas) Q in the point-particle limit, 
solvable at one special value of the dimensionless inverse temperature [3 — 2. 

The situation was changed in very recent years. In a series of works, the bulk and 
surface thermodynamics as well as the large-distance asymptotic behavior of the two- 
body correlation functions were derived exactly for the 2D Coulomb gas of pointlike 
particles, in the whole stability range of inverse temperatures (3 < 2. These results 
were obtained by mapping the Coulomb gas onto the 2D sine-Gordon theory with a 
conformal normalization of the cos-field, and subsequently applying techniques and 
recent achievements in that integrable theory. 

The aim of this short review is to present the exact solution of the 2D Coulomb 
gas to the fluid community, in the language of fluid physics and in a way accessible to 
non-specialists in fleld theory. Unsolved problems and further potential developments 
are pointed out. 
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In section 2, we introduce the 2D Coulomb gas. Section 3 deals with its complete 
bulk thermodynamics, obtained from the mapping onto the bulk sine-Gordon field 
theory. Surface thermodynamic properties of a semi-infinite 2D Coulomb gas in 
contact with an impermeable (ideal dielectric or ideal conductor) wall, and the 
corresponding mapping onto integrable boundary sine-Gordon models are given in 
section 4. The large-distance asymptotic behavior of two-body correlation functions 
in the 2D Coulomb gas is presented in section 5. Section 6 is devoted to miscellaneous 
topics and perspectives. 

2. Basic facts about the 2D Coulomb gas 

The symmetric Coulomb gas, defined in an infinite 2D space of points r G R^, consists 
of point particles {i} of charge {qi — ±1} immersed in a homogeneous medium of 
dielectric constant = 1. The interaction energy of particles is J2i<j QiQjvUri — rj\), 
where the Coulomb potential v is the solution of the 2D Poisson equation 



Explicitly, i;(r) = — ln(|r|/L) where L is a length scale. This definition of the Coulomb 
potential in 2D maintains many generic properties (e.g. sum rules) of "real" 3D 
Coulomb fluids with the interaction potential v{r) — l/|r|, r £ R^. 

The model is treated as the classical one in thermodynamic equilibrium, via the 
grand canonical ensemble characterized by the (dimensionless) inverse temperature (3 
and by the couple of particle fugacities z+{r) = z_(r) = z. We set the free length scale 
L to unity for simplicity; the true dimension of the rescaled z is then [length]'^/^"^. 
The grand partition function is defined by 



where iV+ (N^) is the number of positive (negative) particles and N = + N^. 
Many-particle densities are generated from S in a standard way as functional 
derivatives with respect to the fugacity field Zq{j^), taken at the constant z^{r) — 
Z-{r) — z. At one-particle level, one introduces the number density of particles of one 
sign Uq — Sq,qiH^ ~ ^i)) ■ ^uc to thc chargc symmetry, n+ = n_ = n/2 {n is the 
total density of particles). At two-particle level, one introduces the two-body densities 
'^g<;'(k~i"'l) — iJ2i^j '^9,9i^(''~ ~ I'j))- I* is useful to consider also the pair 
distribution functions gqq'{r) = nqqi{r) / {uqUqi), the (truncated) correlation functions 
hqqi(r) = gqq'{r) — 1 and the Ursell functions Uqq'{r) = nqUqihqqi 

The underlying system of pointlike particles is stable against the collapse of 
positive-negative pairs of charges provided that the corresponding Boltzmann factor 

is integrable at short distances in 2D, i.e. for j3 < 2. To cross the collapse 
point j3 = 2, the pure Coulomb interaction has to be regularized by a short-distance 
repulsion, e.g. a hard-core potential of diameter a around each particle (the particular 
choice of the short-distance regularization influences the results quantitatively, but not 
qualitatively). For small values of the dimensionless density ncr^, the system remains 
in its conducting phase (an external charge is perfectly screened by the system charges) 
up to the Kosterlitz-Thouless (KT) transition of infinite order Q at a specific density- 
dependent Pkt', Pkt = 4 in the low-density limit. In the insulating phase /3 > (3kt, 
the system charges form dipoles and no longer screen an external charge. At high 
enough density, the KT critical line splits into a first order liquid-gas coexistence 



Av{r) 



2nS{r). 



(1) 
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curve 0. In what follows, we shall restrict ourselves to the point-particle Coulomb 
gas in the stability region (3 <2. 

A complete exact analysis can be done in two cases: in the high-temperature 
Debye-Hiickel limit (3^0, and just at the collapse point /3 = 2 Q which corresponds 
to the free-fermion point of an equivalent 2D Thirring model. Although, at a given z, 
the free energy diverges, Ursell functions are finite at /3 = 2. 

As concerns an exact information valid in the whole stability region /3 < 2, 
through a simple scaling argument, the exact equation of state for the pressure P, 
PP = n(l — /3/4), has been known for a long time While the density derivatives 
of the Helmholtz free energy, like the pressure, can all be calculated exactly, the 
temperature derivatives, like the internal energy or the specific heat, are nontrivial 
quantities. Their evaluation can be based on an explicit density-fugacity (n, z) 
relationship. The latter was constructed systematically around the (3^0 point by 
using a bond-renormalized Mayer expansion in density the original bonds —/3v{r) 
are summed up in series to produce the renormalized bonds of strength —f3KQ{Kr), 
where Kq denotes a modified Bessel function and k = (27r/3n)"'^/^ is the inverse Debye 
length. The cluster integrals converge in the renormalized format. The first few 
integrals imply 

^ = 2/^^^^-p I [2C + m (I)] f + ^C(3) (0% C(3) 0(^^)} (3) 

where C is the Eulcr number and ( denotes the Riemann zeta function. 



3. Bulk thermodynamics 



The 2D Coulomb gas is equivalent to the 2D sine-Gordon model ||T^. Introducing 
the microscopic charge density p{r) — X^iLi Qi^i^ ^ r.^), the interaction energy can be 
written as 



E=^J dVdVp(r)w(|r-r'|)p(r') 



-Nv{0). 



(4) 



Let us forget for a while that v{Q) diverges, and rcnormalize the fugacity by the 
self-energy term exp[— /3d(0)/2], without changing the z-notation. Using the fact that 
— A/(27r) is the inverse operator of i;(r) [see equation (|l])], the grand partition function 
(H) with Zq(r) = z can be turned via the Hubbard-Stratonovich transformation into 

/ X>(/)exp(-S'(z)) 



where 



E{z) 



S{z) 



dV 



2z cos{b(f>) 



1/2 



(5 a) 



(56) 



is the Euclidean action of the 2D sine-Gordon theory. Here, 0(r) is a real scalar 
field and / denotes the functional integration over this field. The many-particle 
densities are expressible as the averages over the sine-Gordon action as follows 

riq = zqie'"^"^), nqg,{\r - r'|) = (e'«''*«e''''''*('-')) (6) 

etc. The parameter z in (p6[), i.e. the fugacity renormalized by the diverging self-energy 
term, gets a precise meaning when one fixes the normalization of the coupled cos-field. 
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In the Coulomb system, the behavior of the two-body density for oppositely charged 
particles is dominated at short distance by the Boltzmann factor of the Coulomb 

potential, (r,r') ~ z+Z-\r — r'|~'' as |r — r'| 0. With regard to (^, the 

mapping is supplemented by the short-distance normalization 

^gi60(r)g-ifc0(r)^ ^ |r - r'r*''' as |r - r'| -> (7) 

which was usually omitted in the statmech literature. Under this short-distance 
normalization, the divergent self-energy factor disappears from statistical relations 
calculated within the sine-Gordon representation. This can be easily verified in the 
Debye-Hiickel limit /3 ^ 0, when cos(6(/)) ~ 1 — b'^(fP/2, and the consequent Gaussian 
field theory reproduces the n, z relation (|^) up to the linear /5-term in the exponential. 

In the classical limit of the sine-Gordon theory, only such configurations of the 
(j)-tie\d are considered which fulfill the equation of motion SS = 0. This classical limit 
is integrable i.e. there exists an infinite sequence of conserved quantities. Due to 
the discrete symmetry (p ^ <f> -\- 27rn/& (n integer), the model has an infinite number 
of vacua at </)„ = 27rn/6. The basic "particles", the soliton S and antisoliton S pair of 
equal masses M, interpolate between two neighboring vacua. The S— S pair can create 
bound states, called the breather particles {B}. The sine-Gordon model is integrable 
at the full " quantum level" (all configurations of the c/i- field are considered) as well , 
with the same particle spectrum. The essential difference between the classical and 
quantum theories is that the breathers become quantized, {Bj;j = 1,2,... < 1/^}, 
and their number depends on the inverse of the parameter ^ — 6^/(1 — b^). The mass 
of the i?j -breather is given by 

m, =2Afsin(^^j^ (8) 

and this breather disappears from the spectrum just when mj — 2M. Note that 
breathers exist only in the stability region of the Coulomb gas < b^ < 1/2 
(0 < /3 < 2). The lightest i?i-breather disappears just at b'^ — 1/2 {(3 ~ 2), which is 
the field-theoretical manifestation of the collapse phenomenon. 

Using the Thermodynamic Bethe ansatz, the dimensionless specific grand 
potential 

1 77? ^ 

lim lnS = — ^ (9) 
v^oo V 8sm(7rO 

was found by Destri and de Vega Under the conformal normalization (^, the 

relationship between the soliton mass M and the fugacity z was established in , 

^ n 2-26^ 

(10) 



r(62) 



7rr(l -62) 



2r(^/2) 

where T stands for the Gamma function. Equations (§)"(|lO|) constitute a complete 
set to be solved for the exact n, z relationship 



^/"_.^^/?V^" r(i-/3/4) 
) r(i-h/3/4) 



' 2' 4-/3' 2(4-/3)' 



-/3/4 



(11) 



where F = ^Fx is the hypergeometric function. The expansion of the rhs around 
/3 ^ reproduces correctly the first terms of the renormalized expansion (^). For 
fixed z, the particle density given by ( [III ) exhibits the expected collapse singularity 
n ~ 47rz2/(2— /?) as — > 2^. This behavior can be derived independently by combining 
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an electroneutrality sum rule, — (/riq — X)o'=± l' /"^^^ggK^)? with the short-distance 



asymptotic behavior of n-| (?■) discussed above. Since the derivation of formula ( |10[ ) 

was based on special analyticity assumptions, the check of the results from both sides 
of the stability interval is important. Without noticing it, such checks are made for 
all presented results. 

Based on the explicit n, z relation, one can pass by using the Legendre 
transformation from the grandcanonical to the canonical ensemble, to obtain the 
Helmholtz free energy. We present the explicit result for the excess specific heat 
at constant volume per particle 



0^ 

4 4-13 
2/32 - 



2 r 



16 



V/ 1 - 



(4 - 



f3 



-V/ 

8-/3 
8-2/? 



1 



Att^P^ cos (tt/?/ (4-/3)) 
(4-/3)3 gi^2 



(12) 



Here, ipix) — d[lnr(x)]/dx is the psi fmiction and tp'{x) = X^i^o + *)^- 
/3 — > 2~, the /ks exhibits the expected singularity of type 2/(2 — . Notice that 
the specific heat is independent of the particle density, which is a peculiarity of the 
2D Coulomb gases. 



4. Surface thermodynamics 

Let us now consider a semi-infinite 2D Coulomb gas in the Cartesian half-space a; > 0, 
in contact with a hard wall of dielectric constant in the complementary half-space 
a: < 0. The presence of the dielectric wall manifests itself through particle images ]l^ : 
the particle of charge q at position r = (x, y) has the image of charge q* (dependent 
on cw) at r* = (— x, y). We will consider two particular cases: the ideal dielectric wall 
{ew — 0) with image charges q* = q and the ideal conductor wall {ew — *■ 00) with 
image charges q* = —q. Let us introduce the microscopic charge plus image-charge 
density /5(r) = Qi l^i^ ~ ^i) ^ ^i^ + ^i)] ^{y~yi)\ hereinafter, the upper (-I-) sign 

corresponds to tw = and the lower (— ) sign to ew 00. The interaction energy of 
the particle-image system can be written in both cases as 

E^^Jd^l dV'/}(r)«(|r - r'l)p(r') - ^Nv{0) (13) 

where the integrations over r and r' are taken over the whole 2D space. 

The form of the interaction energy (|l^) resembles the one in (^), and one can 
proceed in close analogy with the bulk mapping. The grand partition function is 
expressible as 

/l?(/)cxp(-g(z)) 
= /P0exp(-5(O)) ^^^'^^ 
where the (/)-field is defined in the whole 2D space and the nonlocal action reads 



S{z) = / d 



l2„ 



^(V0)2 - 2z cos (^-^ [<l>{x, y) ± 0(-a;, y)]^ 



(146) 




To make this field theory local, we introduce two new fields 

[0(a;,2/) +0(-a;,y)] , 0o(a;, y) = [0(a;, y) - 0(-a;, y)] (15) 
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defined only in the positive a: > half-space. The "even" field has a Neumann 
boundary condition, dx<pe\x=o — 0, and the "odd" field a Dirichlet boundary condition, 
<^„U=o = 0. Since J d^r{V^y' = d^r [(V<^e)' + (V(/.o)2] , the fields and 
(j)o are decoupled in the action ( |146[ ). The field, contributing only by its free-field 
part, disappears from ( [14qD by numerator-denominator cancellation. Consequently, 
renaming the kept field by (f>, we arrive at 

/ 2)(^exp(-S'(z)) 



S(z) 



S{z) 



dV 



2z cos(&(/)) 



(16) 



/I?0exp(-5(O))' ^,>o 

Here, the 0- field has the boundary condition: 4)\x=o = for the ideal conductor 
wall; dx4>\x=a — for the ideal dielectric wall. The mapping onto the boundary 
sine-Gordon model is supplemented by the short-distance normalization (uu. Both 
cases under consideration belong to the integrable boundary field theories |l5|]. The 
thermodynamic quantity of interest is the surface tension 7, which characterizes the 
surface part of the grand potential Q = —f3~^ In 2. 

The problem of the ideal conductor wall was solved via a lattice regularization of 
the boundary sine-Gordon model [ p^ , namely the XXZ Heisenberg quantum chain in 
boundary magnetic fields. The surface tension was obtained in terms of the soliton 
mass M as follows 

-1 ^ 



cond 




tan 



77/3 



2(4-/3) 



7r/3 



2(4-/3) 



(17) 



The surface collapse is governed by the interaction Boltzmann factor of a particle with 
its self-image, x~^^'^. The ID integral / Axx'^^^'^ diverges at short distances at point 
(3 — 2 (which coincides with the bulk collapse point), and this is indeed the radius of 
convergence of 7cond- 

The problem of the ideal dielectric wall was solved by exploring a "reflection" 
relationship between the Liouvillc and sine-Gordon theories p7| . The result is 



/37dicl = — 



1 + tan 



7r/3 



2(4-/3) 



cos 



7r/3 



2(4-/3) 



(18) 



At /3 = 2, 7diei keeps a finite value |l8j. The analytic continuation of ( |18| ) beyond the 
bulk collapse point predicts a surface collapse at /3 = 3. Such a phenomenon is due to 
a paradoxical short-distance attraction of a particle with its own image charge of the 
same sign [|l^ . 

The surface thermodynamics of a plain hard wall {ew = 1), which is the last 
exactly solvable case at the free-fermion point [3 — 2 ||l^, is an open problem. 



5. Large-distance behavior of particle correlations 

In a 2D integrable theory characterized by a "particle" spectrum, correlation functions 
of local fields can be written as an infinite convergent series over multi-particle 
intermediate states, in terms of the corresponding form-factors. The form-factor 
representation is especially useful at large distances, since the dominant contribution 
of the series comes from an intermediate state with the minimum value of the total 
particle mass. In the sine- Gordon theory, for topological reasons, solitons S and 
antisolitons S coexist in pairs, the total mass of the pair being 2M . The breathers 
{Bj}, when they exist, have lighter masses [see formula (||)]. The lightest Bi- 
breather with mass mi = 2Msin(7r^/2), which exists in the whole stability region 
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of the Coulomb gas, governs the large-distance asymptotic behavior of the two-point 
correlation in (^. In particular, one has for the correlation function h^qi [ |20[ |, as 
r oo, 

hqq'i"^) ^ m'KT), ^^''''"^^(,2^^) cxp(-mir) (19) 

where A is a /3-dependent prefactor. The specific dependence of hqqi (r) on the charge 
product qq' means that the two-particle correlations are determined at large distance 
by the charge-charge correlation function hp{r) = (1/4) ^,^_|_ qq'hqq'{r) = h{r). 

On the other hand, the density correlation function /i„ = (1/4) g'=± ^99' 
vanishes for the lowest one-i?i breather state, and becomes nonzero only for the next 
two--Bi breathers state The mass of two i3i-breathers, 2mi, is smaller than the 
one of the soliton-antisoliton pair, 2M, in the region < /3 < 1. Consequently, as 

, -.-V - . . forO</3<l ; 
hJr) (x{ , \ 20 
«V J ^ , for 1 < /3 < 2 . ^ ' 

The correlation length depends continuously on /?, but its first derivative with respect 
to /3 is discontinuous at /3 = 1. The large-distance exponential decay of /i„ is faster 
than the one of hp. The two correlation lengths coincide just at the collapse point 
/3 = 2, where hp and hn differ from one another only by the inverse-power law 
prefactors. 

6. Miscellaneous topics and perspectives 

The 2D Coulomb gas exhibits, at any temperature in the stability region, a universal 
finite-size correction to the grand potential as if we had a critical theory with the 
conformal anomaly number c — —1, although the particle correlation functions 
presented here decay exponentially. This phenomenon follows intuitively from the 
sine-Gordon representation of S (p^, with the critical massless Gaussian field theory 
(c = 1) in the denominator. The explicit checks of the critical- like behavior were done 




at the free-fermion f3 — 2 point for various geometries of confining domains |22, 23|, |24| , 
and at any /? < 2 for the sphere ||2^, |2^ and for the disk . 

The ultimate task is to solve exactly the 2D Coulomb gas with a short-distance 
(maybe temperature-dependent) regularization of the pure Coulomb potential. We 
have made a first step toward this aim by deriving the leading correction to the exact 
bulk thermodynamics of pointlike charges due to presence of the hard core of diameter 
a around particles |^. The results, which are conjectured to be exact in the low- 
density limit up to /3 = 3, reproduce correctly the cr-singularities of thermodynamic 
quantities at /3 = 2. They also confirm a "subtraction" mechanism of singularities 
between the collapse point p = 2 and the KT transition point /3kt — 4 within an 
ansatz proposed by Fisher et al. |28| (excluding the existence of an infinite number of 
intermediate phases proposed in ||29|), however, predict a different analytic structure 
of this ansatz. There exist candidates among integrable ID quantum systems, for 
example the lattice sine-Gordon model which, after being formulated as a 2D 
Euclidean theory, might represent a Coulomb gas regularized in the whole temperature 
range. 

Another topic, which attracted much attention in the last years due to the 
phenomenon of charge inversion plj], is the charge asymmetry. The Coulomb gas 
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with \q^\ — 2\q_\ was solved in terms of the equivalent complex BuUough-Dodd field 
theory in reference |Q. Here, the fundamental changes in statistics caused by the 
charge asymmetry (e.g. the same correlation length for both correlation functions hp 
and hn at any stable /3) were documented. This result might be a motivation for 
the exact solution of the 2D one-component plasma, which is the extreme charge- 
asymmetry case of the 2D Coulomb gas. 
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